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Abstract. Universal limits for the eigenvalue correlation func- 
tions in the bulk of the spectrum are shown for a class of non- 
determinantal random matrices known as the fixed trace or the 
Hilbcrt-Schmidt ensemble. These universal limits have been proved 
before for determinantal Hcrmitian matrix ensembles and for some 
special classes of the Wigner random matrices. 

1. Introduction and the statement of the result 

Let 7ijs[ be the set of all iVx N (complex) Hermitian matrices, and let 
txA = Ylii=i a a denote the trace of a square matrix A = (aij)^j =1 . Hn 
is a real Hilbert space of dimension N 2 with respect to the symmetric 
bilinear form (A,B) i— ► tiAB. Let denote the unique Lebesgue 
measure on Ti N which satisfies the relation In(Q) = 1 for every cube 
Q C T~Cn with edges of length 1. A Gaussian probability measure on 
Ti.N invariant with respect to all orthogonal linear transformations of 
H.N is uniquely defined up to a scaling transformation. Such measures 
form a one-parameter family (fi s N ) s >o, where the measure fj, s N is specified 
by its density 



with respect to 1^. Thus, for a random matrix X distributed according 
to n s N we have 



The set Hn endowed with the measure /j, s n is called the Gaussian Uni- 
tary Ensemble (GUE). Let X be a random N x N Hermitian matrix 
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(that is a random variable taking values in TCn). We consider the eigen- 
values Ai, A 2 , . . . , Ajv of the random matrix X as a finite sequence of ex- 
changeable random variables. By definition, this means that their joint 
distribution P§ does not change under any permutation of these vari- 
ables. Let for each n, 1 < n < N, P^ N denote the joint distribution of 
some n of these N variables. Obviously, P* N is a permutation invariant 
probability measure in IR n . In particular, the measure Pf N describes 
the distribution of a single eigenvalue. By definition, the n-point corre- 
lation measure p* N of a random matrix X is a non-normalized measure 
defined by the relation 

( L3 ) Pn,N = ( N _ n y P n,N- 

For a measurable set id" the quantity p^ N (A) can be interpreted 
as the average number of n— tuples of eigenvalues in the set A. If the 
measure p^ N is absolutely continuous with respect to the Lebesgue 
measure on IR n , its Radon-Nikodym derivative R*n * s called the n- 
point correlation function of the random matrix X. In particular, the 
measure p* N has total mass N. For a measurable set E C IR 1 , the 
quantity pf N (E) expresses the expected number of the eigenvalues be- 
longing to E. The corresponding density with respect to the Lebesgue 
measure in M 1 , if it exists, is called the eigenvalue density or the density 
of states (caution: under the same names the normalized versions of 
the same measures are considered in the literature as well). Let X N be 
a random matrix with the distribution p s N . For n — 1, . . . , N we set 

pGUE,s = p Xx and p GUE,s = p X^ A classical regult for the GIJE says 

that we have 

(1.4) P?""- W, 

where the measures converge in the weak sense, and W is the standard 
Wigner measure on [—2, 2] defined by the density 



(1.5) w(x) = (2tt)-V( 4 -^ 2 )+> xeR. 

In terms of the correlation measures the same relation reads 

(L6) N P ^ n^ W - 

For the n— point correlation measures we have a similar relation 

(1-7) jL-Sr'V 

iV" ' N^oo n times 

which means that the eigenvalues become independent in the limit. 
However, for n > 2, the study of a finer asymptotics near a point from 
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the principal diagonal in the cube (— 2,2) ra shows [TUj [TT] : for every 
u e (-2,2) and h, . . . ,t n G K 1 



;i.s) 

1 



ji^ , Nw , u ^ n RnT N ( u + {h/Nw(u)), ...,u+ (t n /Nw(u))) 



det 



/ sin 7r(ij — t 



i.j=l 



This limit relation presents a pattern for many other results, in partic- 
ular, for that of the present paper. The right hand side of this relation 
represents an example of the correlation function of a so-called determi- 
nantal (or fermionic) random point process [15]. In general the n— point 
correlation function R n of such a process is given by the formula 

R{ui, ...,u n )= det\K(ui,Uj)\l j=1 , 

where K is the kernel of an integral operator on the line, which is trace- 
class having been restricted to any finite interval in R, and subject to 
some further conditions (see [UJ for a detailed exposition). Moreover, 
in the asymptotic Hermitian random matrix theory (Kn)n>o are the 
reproducing kernels of the subspaces of polynomials of degree < N — 1 
with respect to some weight on the line. In this case we call the corre- 
sponding matrix ensemble determinantal. The GUE gives an example 
of such an ensemble. To a large extent the asymptotic study of deter- 
minantal ensembles reduces to that of the respective kernels [2J [5] . I n 
particular, the same local limit as in the case of GUE is established 
in [T2J E] for two broad classes of determinantal matrix ensembles. It 
is expected that the sin-kernel limit (first discovered by F. Dyson) is 
rather common. This is known as the universality conjecture. Outside 
the class of determinantal Hermitian random matrices only very few 
results on the asymptotics of the correlation functions are known (see, 
for instance, [10], where a mixture of determinantal measures is con- 
sidered) . 

In the present paper we investigate the following non-determinantal 
ensemble of Hermitian random matrices. Let 

(1.9) S r N = {A e H N : tr A 2 = r 2 } 

be the sphere in Hn of the radius r > centered at the origin. Set 
r = ^N. The sphere S^ N carries a unique probability measure v s N 
invariant with respect to all orthogonal linear transformations in the 
space T~Cn- We call this measure the fixed Hilbert- Schmidt norm ensem- 
ble (or just HSE) to reserve the term "the fixed trace ensemble" for more 
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general ones (see [T]). Let Y N be a random matrix distributed accord- 
ing to i*. We set for n = 1, . . . ,N P^ s = P^ N and p*f s = p^ N . 
It is a well known result (see [TTJ [H]) that 

(1.10) pBSEA/N ^ ^ 

N^oo 

like in the case of GUE. In this paper we prove that the correlation 
functions R^^' 1 ^ of arbitrary order n (1 < n < AT) near every 
point u G (—2,2), have the same determinantal limit with the kernel 
sin7r(ti — t2)/^{t\ —£2) as the GUE correlation functions (for n = 1 the 
limit equals 1). More precisely, we establish in this paper the following 
result. 

Theorem. Let R^n ^ e ^ e n ~ P ^ correlation function of the eigen- 
values for a random matrix uniformly distributed on the sphere S^. 
Then for every u G (—2, 2) and ti, t 2 , ■ ■ ■ , t n G M 1 
(1.11) 

as iV — > 00. For even/ a G (0,1) and A > £/ie relation (CQ) ZioZds 
uniformly in all u G [—2 + a, 2 — a], and ti G [—A, A], . . . , t n G [—A, A]. 

The class of fixed trace matrix ensembles was studied in several re- 
cent publications. The authors are indebted to G. Akemann for draw- 
ing their attention, after the first version of the present work has ap- 
peared as a preprint, to the paper [2J, he published jointly with G. 
Vernizzi. In this paper the local universality is studied for a class of 
ensembles containing the HSE. The authors provide heuristic argu- 
ments for the universality near zero based on the complex inversion 
of the Laplace transform. Our proof is based on the complex Laplace 
inversion as well (more precisely, we arrive at the Fourier transform 
after a sequence of translations of the origin and rescalings depending 
on N). However, establishing bounds for orthogonal functions and re- 
lated kernels in the complex domain allows us to rigorously conclude 
convergence near every point of the interval (—2,2). 

In the following we will sketch the main steps of the proof. The 
guiding principle is that results for the Hilbert-Schmidt ensembles are 
deducible from the corresponding results for the GUE using the 'equiv- 
alence of ensembles' or the 'concentration phenomenon'. In our setup 
a primitive form of the concentration is given by the law of large num- 
bers for the squares of the Hilbert-Schmidt norms of the GUE random 
matrices. Supplemented by some estimates of the probabilities of large 
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deviations, this is the main tool in [9], where the universality is shown 
near zero for the correlation measures (rather than the correlation 
functions) of the Hilbert-Schmidt ensemble; convergence was shown 
there in the weak topology determined by the continuous compactly 
supported functions. Such simple arguments seem to be insufficient for 
proving local results for correlation functions of eigenvalues (maybe, 
with exception for eigenvalues near zero), particularly, for stronger 
topologies. This agrees with M.L.Mehta's doubts ([II], Sect. 27.1, 
p. 490) concerning the deducibility of the local results for correlation 
functions of the Hilbert-Schmidt ensembles from the corresponding re- 
sults for GUE by using the equivalence of ensembles. Solving this open 
problem in the present paper, we use a local form of concentration 
given by the local central limit theorem for the densities of the squared 
Hilbert-Schmidt norms. First we represent the Hilbert-Schmidt mea- 
sure as a conditional measure of the GUE, given the Hilbert-Schmidt 
norm of the GUE random matrix. Starting with the disintegration 
of the GUE according to the level sets of the Hilbert-Schmidt norm, 
we arrive at formula (I2.9P which is the crucial ingredient of the proof. 
Here we have to extend the scaling parameter to the complex domain. 
The formal Fourier inversion applied to this formula gives an 'heuristic 
proof of the result. For an outline of it see Section 2. To make this 
sketch rigorous we need asymptotic estimates in the complex domain 
for the kernels related to the Hermite functions. This is done in Section 
3, based on the results in [3], [I], [5]. Unfortunately, some of the results 
we need are contained in these papers not explicitly enough and have 
to be extracted from the proofs rather than from the statements (see 
the proof of Lemma ED- The main convergence result we need is Lemma 
121 All necessary bounds are summarized in Proposition [U Based on 
these two facts, we complete the proof. Note that the analytic part of 
the present paper may be viewed as a form of the Tauberian theorem. 

The authors are indebted to the anonymous referees for careful read- 
ing of the manuscript. Moreover, the main theorem in the first version 
of the paper required excluding of a neighborhood of zero (this case has 
been considered in the note [9]). Stimulating questions of one of the 
referees and the editor led us to a revision of the proof which removes 
this restriction. Finally, we would like to thank Justine Swierkot for 
her assistance in the preparation of the text. 
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2. Disintegration, a Fourier transform formula and the 

sketch of the proof 



In this section we discuss a disintegration representation of the GUE 
in terms of the HSE, and derive a Fourier transform formula involving 
these matrix ensembles. We suppress in this section the "spectral" 
arguments of the correlation functions and related quantities assuming 
that these arguments vary inside the domain described in Section 1. 

For every r > denote by S N the sphere of radius r inKjv centered 
at the origin. Let for s > be a random matrix in Ti^ distributed 
according to /% UE ' S . Set T/v = tr(X N /s) and Yn = NXn/ y/trX^. 
The random variable T/v can be represented as a sum of N 2 squares 
of independent standard Gaussian random variables, hence it follows 
the familiar \n 2 distribution. Moreover, T N and Y N are independent, 
and Y N is uniformly distributed on the sphere in TCn, that is, Y N 
is distributed according to v N in the notation of the previous section. 
Then can be represented as 

Yj 



(2.1) X N = ^^T N ~ 

with T/v and Yn as above. Let 7^2 denote the the probability density 
of T/y. Then it follows from (12.11) that 

roc 

(2.2) = u n N lN-*{s~ l u)s~ l du. 

Jo 

As a consequence of (12.21) . the correlation functions of GUE and HSE 
for 1 < n < N — 1 satisfy the relation 

(n o\ D GUE,s / t-jHSE,m/JV i -\ \ -lr 

I 2 - 3 ) R n.N = / R n,N 1n 2 [ s U)S du. 

JO 

Note that 



ET N = N\ DT' N = E(T N - ET N ) 2 = 2N 2 
and, for every m > 0, 

' {2 m / 2 T{m/2))- 1 x^- 1 e- x/2 , ifu > 0, 



(2.4) lm {u) 



0,if« < 0. 



Set for s > 7 m , s (-) = s _1 7 m (s -1 -), so that 7 TO>1 = 7 m (note that 
the same set of densities with a different parametrization appears in 
Lemma [8] as (f a , P )a, P >o)- Observe now that the probability density of 
sTjv is given by 7jv 2 «(•)■ Then (12.31) can be rewritten as 



00 

2 



(n r-\ D GUE,s / D HSE,u/N- , s , 

(2.5) R„ N = / R N -f N 2 tS (u)du. 







In particular, we have 

fn c x D GUE,1/7V / D HSE,u/V 2 / x , 

I 2 - 6 ) R n,N = / R n.N lNM/N{U)du. 

JO 

Our goal is to investigate the limiting behavior of R^^' 1 ^ when n is 

fixed, N — > oo, and the "spectral" arguments of -R^^' 1 ^ vary within 
an e/iV— neighborhood of a point from the main diagonal in the cube 
(—2, 2) n . However, we prefer to consider the function u \— > R^ S ^' U ^ N 

rather than its value R^jy' 1 ^ at AT. Moreover, we will perform the 
study of the latter function indirectly, first dealing with the product 
u i — ► R^ S ^' U ^ N 7n 2 ,i/n( u )- After the change of variable (u — N)/y/2 = v 
in (12.61) we obtain the relation 

/oo 
-oo 

where 

(2.8) g^) = R^ 1/N+vV ~ 2/N2 V2 lN ,, 1/N (N + vy/2). 

Notice that v i— > 7n 2 ,i/n(N + vy/2) is the probability density of 
the centered and normalized random variable (T^a — N 2 )/V2N 2 , and 
it tends to the standard normal density tp : v i— ► (l/\/27r) exp (— f 2 /2) 
as iV — ► oo. Thus, the limit behavior of the density 7n 2 ,i/n{N + ■ y/2) 
is well understood, and we have to study q^i')- For every fixed u, 
due to the relation 7^2 ,,(•) = s -1 7jv2(s -1 -), we can analytically extend 
7jv2 jS (w) to the the domain 9fts > 0. Moreover, in the formula (12.51) 
the integral in the right hand side can be analytically continued in s 
in accordance with the continuation of "frf* a (u) mentioned above. This 
leads to the corresponding continuation of i?^ E,s so that (12.51) holds for 
s from the right half-plane. Now we will evaluate the Fourier transform 
of the (nonprobabilistic) density q N 2 keeping in mind that by (12. 7p and 
(12.81) it is a nonnegative integrable function. 

Lemma 1. 

POO 

(2.9) / exp {ipv)q N *{v)dv = ^(p)R^ W -^ /N)N) , 



where 

(2.10) <j) N 2(p) = exp (-ipN/ V2)(l - ip^2/N)-^ 2/2) 



is the characteristic function of the random variable (T^2 — N 2 )/\ / 2N 2 . 
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Proof. Denoting by C m the normalizing constant in formula (12.41) . we 
have 
f2.ll) 

o 

exp {ipv)qjq2{v)dv 

K> 
O 

exp ( IP v)R^ 1/N+vV ~ 2/N V2 lN2A/N (N + vV2)dv 



. ip(u -N)\ HSE, U /JV 2 , . , 

exp — )R n>N 'j N 2 A/N (u)du 



V2 

poo 

exp (^F) / exp (^ u "^2)- R ^'" /iV lNM/N(u)du 

ex P (if) r Rnf^C^NiNu)^- 1 exp (-^ (l - ^))d« 



V+ 1 



-p(^)(i-¥) "■ >< 



N z -. 

J^^C^N (Nu (l - ^)) 2 "exp (-^ (l - ^f))du 



ex P v ^r) ( X - iJ #) 2 / w <^ U/ ^7 wa ,i/( ( i-W5/ivw(«)^ 
^iy 2 (p)<S E ' 1/((1 -^ /iV)JV) . 

□ 

In the following we will outline our approach. Write 

1 , n s 1 f°° , / x 1 D GUE,l/((l-ip v / 2/A r )A r ) j 

:Qn^{0) = — <p N 2{p) u A jv d P- 



(Nw(u)) n ^ v ' 27ry_ 00 rJV ^'(JViufa))' 

Passing to the limit in the integral on the right hand side (uniformly 
with respect to the spectral variables), the right hand side has the same 
limit as 

or, in view of the local Central Limit Theorem (CLT) , as 

1 pGUE,l/iV 

spM{Nw{u)) n n ' N 
On the other hand, it follows from (12 .81) that 

(2 - 12) jN^r qmi0) = W^^ 1,NV ~ 2 ^ {N) - 
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Again, the local CLT implies that \f2^ N 2^/ N (N) - ► \ j\phx. Therefore 



1 



-R 



HSE,1/JV 



(Nw(u)) n n ' N 
tends to the same limit as 



pGUE,l/AT 



(Nw(u)y 

and the conclusion follows. In the next Section these heuristic argu- 
ments will be made rigorous. 

3. Proofs 

For every a G R, through the rest of the paper, we will denote by 
(■) a the function 

(3.1) (■)<* : C\ (-oo,0] -> C : z i-> exp a log z, 

where log denotes the principal branch of the logarithm. We use as 
a notation for (-) 1 / 2 extended to by v^O = 0. 

The same symbol (for example, C, L, M, R, a, ... ,5 with or without 
indices, Nq, and so on) may denote different constants in the bounds 
obtained in the rest of the paper (even in the same proof). However, we 
will indicate explicitly the dependence of such constants on parameters. 

First we are going to state some definitions and basic formulas re- 
lated to the Hermite polynomials and Hermite functions. Then we will 
formulate some results on the asymptotics of the Hermite polynomials 
and Hermite functions in the complex plane. The asymptotic behavior 
of the Hermite polynomials at the scale we are interested in was first 
established in 1922 by Plancherel and Rotach. For a convenient form 
of these (and much more general) results we refer to the monograph 
[5] and the papers [3J,[4] (in particular, Appendix B), and j6]. These 
results will be employed later in this section. 

Let s > 0, and let for every N, N > 0, Pn(', s ) be a polynomial of 
degree N with a positive leading coefficient such that for (Pn(-, s))n>o 
the relations 



Pm(x, s)Pn(x, s) exp (— x 2 / (2s)) dx = Sm,n, M, N = 0, 1 



hold. 

The (monic) Hermite polynomials (Hn{ 4 , s))jv=o,i,... are defined by the 
expansion 

(3.2) exp ( 7 o; - 7 V2) = £ H N (x, s)^- 

N=0 
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and for M, N = 0, 1, . . . satisfy the relations 

H M (x, s)H N (x, s) exp (-x 2 /(2s)) dx = 8 m ,n^s^ m+n+1)/2 N\ 
so that 

M-> *) = (27r) l/4 g (27V 1 + l)/4 (iV!) l/2 ^(-' S), JV = 0, 1, . . . 

The polynomials Ptv(-,s) satisfy the difference equations 
(3.3) 

xp N (x, s) = s 1/2 VN + lp N+1 (x, s) + s 1/2 p N -i{x, s), N = 1, 2, . . . 

The Hermite functions (<Pn(', s ))n=o,i,... defined by 

fi N (x, s) = p N (x, s) exp (-x 2 / (4s)), N = 0, 1, . . . , 

form an orthonormal sequence of functions in L 2 (K, A) where A is 
the Lebesgue measure in KL Let us introduce the standardized Her- 
mite polynomials and Hermite functions (corresponding to the weight 
exp (— x 2 /2)) by the relations 

H N (-) = H N (; 1), p N (-) = p N (; 1), tp N (-) = tp N (; 1) 

so that 

Hn(; s) = s N / 2 H N {- s- 1 ' 2 ), p N (; s) = s-^p N (. s^ 2 ), 

Vn(;s) = s- 1 ^n(-s-V 2 ). 

The above relations for the Hermite polynomials imply that the Her- 
mite functions satisfy, for every k > 1, the following system of differ- 
ential equations: 

¥k( x ) = ~n ( Pk( x ) + Vk(p h -i(x), 

( 3 - 5 ) , r x 

= -Vkip k (x) + -(fk^ix). 

The reproducing kernel K N {-, -, s) of the orthogonal projection in L 2 (1R, A) 
onto the linear span of 0o(-, s), . . . , <£>jv-i( - > s ) is given by 

JV-l 

(3.6) K N (x, y,s) = ^2 <Pk(x, s)<p k (y, s). 

k=0 

Note that the correlation functions of the GUE can be expressed in 
terms of the latter kernel as 



(3.7) R®™ ,8 (xi, ...,x n ) = det (k N (xi, x j} s)j 

(see PU, (6.2.7)). 
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Further, setting 

K N (x,y) = K N (x } y } 1), 

we obtain 

(3.8) K N (x, y, s) = s- 1/2 K N {xs- ll \ y S - l/2 ). 

The Christoffel - Darboux formula for the kernels Kn and Kn reads 

f o n \ v- ( \ 'PN(x,s)<p N - 1 (y,s)-<p N - 1 (x i 8)0 N (y,s) 

(6.9) K N {x, y, s) = 

x-y 

and 

(3.10) K N (x,y) = . 

x-y 

Note that the kernels Kn, Kn enjoy the property 

(3.11) K N (-zi, -z 2 , s) = K N (z 1 , z 2 , s), K N {-z u -z 2 ) = K N (z u z 2 ). 

As to K N , this is a consequence of (I3.1(jp and the fact that the functions 
(<£>jv(-))iV>o are even or odd depending on N being even or odd; for 
Kn it follows from (13. 8p . The following integral representation of the 
reproducing kernel is a version of formula (4.56) in [8]: 

K N (x,y) = 

(3.12) 



l~N f°° 



(y + r) + (f N -i(x + r)(p N (y + r))dr 



so that 
(3.13) 

K N (x,y,s) = 

IW f°°/ \ 

Y s I K(*+^)^-i(^+r)+^-i(^ + r)^(^ + r))dr. 

The relations (13.41) . (13.81) . and our extension of the function (■) a to C \ 
(— oo, 0) allow us to continue H N (x, •), Pn(%, •)> ^(x, ') an d K N (x, y, •) 
to this domain analytically in the parameter s. The relations (13.21) - 
(13. lip remain valid under these continuations. So does (13 . 13[) whenever 
the integrals there are well-defined. 

Denote by w the analytic continuation of the standard Wigner den- 
sity (11.51) to the domain C \ ((— oo, 2] U [2, oo)) . Let us define for every 
a G (0, 2) and (3 > the sets S a> p and S a ^ by the relations 

S a>/3 = {zeC: \Uz\ <2-a, \Qu\ < (3} 

and 

S^p = {zEC: \$lz\ < 2 - a, \%u\ < f3}. 
n 



Set for H G R d(H) = yl+iH and observe that 
(3.14) \d(H)\ = {l + H 2 f/\ 



(3.15) Ud(H) 



(3. 16) 3fd(jy) = sgn # ' ' 1 
and 

(3.17) ($ld(H)) 2 - (%d(H)) 2 = 1. 
It is clear from (13.171) that 

(3.18) \&(ud{H))\ > \$f(vd(H))\ 

for every uei Note that for every b > the equation 

\&d(H)\ = b 
has a unique nonnegative solution 

(3.19) H b = v/(l + 26 2 ) 2 -l. 

Lemma 2. Let a, /?, and u be real numbers such that a G (0, 1), 
< (3 < a, and u G (—2 + 2a, 2 — 2a). T/ien the following assertions 
hold true: 

(1) For every bE [—00,00] the relations \H\ < \Hf,\ and \Q(d(H)) \ < 
b are equivalent (we set here = 00, 3(<i(oo)) = 00, and 
^s(d(—oo)) = —00); the relations \H\ < \Hb\ and \$s(d(H))\ < b 
are also equivalent. 

(2) For every real H we have ud(H) G S a p whenever \Q(ud(H ))[ < 

(3) If for a certain real H the relation \H\ < \Hp/ 2 \ (or, equiva- 
lent^, \^s[d{H)j\ < [3/2) holds, it follows that ud(H) G S a> p. 

Proof. It is clear from (I3.19P and the definition of H\, that for b G 
[0,oo] and H G [0, 00] the functions b ^ and H \— > are 
strictly increasing mutually inverse mappings. This shows equivalence 
of the conditions in point 1 of the lemma for the functions in question 
restricted to [0, 00]. The same is true for b,H G [—00,00] since both 
functions are even. This proves assertion 1. 
For every H such that 

(3.20) \%(ud(H))\ < 

12 



we have 

(3.21) \&(vd(H))\ < 2 - a 

since it follows from the relation (13. 17ft and the assumptions (3 < \fa, 
u G (—2 + 2a, 2 - 2a) that 



\M(ud(H)) \ = \/u 2 + ($fu(d{H))) < v / (2-2a) 2 + p 2 

< ^(2-2a) 2 + a 2 < 2 - a. 



In view of the definition of S a) p, the second assertion is proved. The 
third assertion is an immediate consequence of the second one since 
|u| < 2. 

□ 

Corollary 1. Let a, f3, andu satisfy the assumptions of Lemma\^ and 
assume that the relation ud(H) ^ S a ^ holds for some real H. Then we 
have \$t(ud(H))\ > (3. Moreover, the relation \50t(ud(H , ))\ > (3 holds for 
every real H' such that \H'\ > \H\. 

Proof. The corollary follows immediately from the second and the first 
assertions of Lemma [2] combined with the fact that, by (13.161) . |9fd(.H")| 
is a strictly increasing function of 

□ 

In the next Lemma the main result on convergence is formulated. 

Lemma 3. For every number a G (0, 1) there exists a positive number 
H = H(a) such that for every A > the relation 

(3 Nh k " ( (- + 14>) dm (- + dm ^) 

[ ' ' sm(ir(ti-t 2 )d(H)w(ud(H))/w{u)) 
n^o n (ti -t 2 )d(H) 

holds uniformly for all real numbers u, ti,t 2 , and H provided that u G 
[-2 + 2a, 2 -2a], [* x | < A, \t 2 \ < A, and \H\ <H. 

Proof. Note that for every a G (0, 1) there exists (3' = (3'{a) > such 
that for every A > the relation 

1 ~ / zi z 2 1 \ sin7r(zi — z 2 ) 

(3.23) ——K N [v + -rr^,v+ ' ] 



Nw(v) V Nw(v)' Nw(v)' N ) N^oo n(zi — z 2 ) 



holds uniformly for all complex numbers v G S a ^t( a ), Z\, and z 2 such 
that \z x \ < A, \z 2 \ < A. 

This assertion (and more general ones involving some class of weights) 
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is implicitly contained in the papers [3] and [I] (see also the mono- 
graph [5j). Actually, Lemma 6.1 in [3] establishes the desired result 
for real u, Z\ and z 2 . However, the same reasoning applies to the com- 
plex u, z±, and z 2 satisfying the assumptions just made provided that 
f3' is chosen to be sufficiently small in accordance with a. The key 
ingredient of the proof in these references, the boundedness property 
of a certain derivative, is established in [3J and [I] for some complex 
neighborhood of an arbitrary real point u G [—2 + 2a, 2 — 2a] (see 
relation (4.122) in [3]) which allows to bound this derivative in a rect- 
angular strip S aj 0't a ). Lowering (3' if necessary, we can substitute f3'{a) 
by some /3(a) < a/A. Note that the function w(-) has no zeroes in 
C \ ((— oo,2] U [2, oo)) D S a> p( a y Picking a point v' G S a ^ a ) and 
setting in (I3.23P z\ = z^wiy) /w{v') (i = 1,2), we obtain 
(3.24) 



-K N \ r : _ Z \ - ,v+ 



1 \ sm I 



Nw(v') V Nw(v')' Nw(v>)' N J n->oo n(z[ - z' 2 ) 

Note that for v G S a ^ a \ we have < |w(f)| < C with some C > 1. 
This fact and (I3.23j) imply that for every A > (I3.24p holds uniformly 
in v,v' G Sa^a) and |^| < A, \z' 2 \ < A. 

According to assertion 3 of Lemma [21 for every u G [—2 + 2a, 2 — 2a] 
and every H G [— H^ a y 2 , Hp( a y 2 ] we have ud(H) G S a ,p( a )- More, it 

is clear from ( 13.141) that for H G [— Hp( a y 2 , Hp( a y 2 ] we have C' < 
|rf(i^)| < C with some positive constant C > 1 depending on a. Then, 
setting v' = u E [—2 + 2a, 2 — 2a], t> = ud(H) G S a) /3( a \, z\ = tid(H) 
(i = 1,2), we obtain @£22). 

□ 

Remark 1 . It is sometimes convenient to use the following reformulation 
of (I3.23P and the related assumptions : 

For every number a G (0, 1) there exists a real numbers (3 = /3(a) < a/A 
such that for every A > the relation 



1 ~ sin(AT7r(wi -w 2 )w(w)) 

:K N {V 1 ,V 2 ,N )- 



— > 

AT— >oo 



AT JVV A ' z ' ' Nir(vi — v 2 ) 

holds uniformly for all complex numbers v, Vi,v 2 satisfying the con- 
straints V G S a> /3( a ), \v\ — v\ < A/N, \v 2 — v\ < A/N. 

We will now derive upper bounds for the modulus of the kernel 

1 

w 



1 :K N (;;N- 1 ) 
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valid for various values of its arguments. We will use asymptotics of 
two types for the Hermite functions: one of them is valid in a nar- 
row strip around the interval [—2 + a, 2 — a]; another one applies in 
the complement to some neighborhood of [—2, +2]. With this asymp- 
totics, after obtaining certain intermediate bounds for the quantity 
\N- 1 K N (u + t 1 /N)d(H),u + t 2 /N)d(H),N- 1 )\ , we will finally deter- 
mine for this quantity an upper estimate which does not depend on 
H G R. This bound and analogous bounds for the correlation functions 
will be derived in Proposition [TJ 

In the next two lemmas we present upper estimates for the Hermite 
functions. These estimates are immediate consequences of the Plancherel- 
Rotakh-type asymptotics. First we consider a bound valid in the strip 
Sa,/3(a) where [3(a) was defined in the proof of Lemma El 

Lemma 4. For any a G (0,1) there exist a constant C(a) > such 
that we have 



for every N G N and z G S a M a y 

Proof. First we recall a known result about the Plancherel-Rotach 
asymptotics for the Hermite functions in a rectangular strip (Theorem 
2.2, part (ii), in [I]; the particular case of the Hermite functions for the 
measure exp (— x 2 )dx is considered there in Appendix B). This strip 
in the complex plane contains, up to some neighborhoods of the ends, 
the interval where the zeroes of the Hermite functions asymptotically 
concentrate. Set 



(3.25) 




and 



(3.26) 




if>:C\ ((-oo, -1] U [1, oo)) -»• C : z i-> 




The function arcsin is defined as the inverse function of 
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In our notation, the statement in |f] reads: for every a G (0, 1) uni- 
formly in z G {z' G C : |9£(V)I < 1 — a, \Q(z')\ < a} we have 

<p N (V2Nz,l/2) 

= (1 - ^)" 1/4 (1 + z)- 1 ^ 



(3.27) 



7P 



'2N 



x jcos^iW ^ ip(y)dy + ^ arcsin ^1 + f^ll 

+ sin^A^7r J i[)(y)dy — - arcsin z\ O (— ^ \. 

The analytic functions z i— ► f*ip(y)dy and arcsin are bounded on 
the set {V G C : < 1 — a, |3(z')| < a}. Representing cos and 

sin as the sum and the difference of the exponentials, it follows from 
(13.27!) (omitting the multiplier iV _1//4 ) that the modulus of the left hand 
side of (13.271) is bounded above by the N— th power of a positive number 
depending on a. For every a G (0, 2), due to the relations 



<p N {zVN) = 2- 1 / V ( -^VN, 1/2 ) = 2- 1 / V ( fv^iV, 1/2 



V2 ' ' ) ™\2 

a similar bound by the N— th power of some number depending on a 
holds for (pisr(z\/ r N) uniformly for z G S aa . For 



ip N ^{zVN) = (pK.^zy/iN-lWN/iN-lj) 

we also have a bound of such type in S , Q , iQ / 4 C S a /2, Q /2- Indeed, set 
Pn = \/ N/ (N — 1), and let N = N (a) be defined by the equation 
Pn = (2 — a/2)/ (2 — a). Then we have pn < 2 and pNS a>a /4 C S Q / 2)a /2 
whenever iVeN satisfies N > N . This means that for z G S a ^ a u both 

^ G S a / 2) a/2 and z^N/(N - 1) G S a / 2 , a /2 hold for every iV > A^ . 
Therefore, the conclusion of the Lemma is valid for such values of 
N. Choosing a larger constant C(a) if necessary and recalling that 
(3(a) < a/4, the proof is completed. 

□ 

We will obtain now some estimates valid the Hermite functions in 
the domain \5sz\ > 5 > 0. 

Lemma 5. For every 5 > there exist constants C{8) and M(S) such 
that the inequalities 



(3.28) \p N (zVN) I < C(8)N~ 1/4 M N (8) \z 
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A' 



and 

(3.29) |^_i(*VjV)| < ^(^(Ar)- 1 /^- 1 ^) l^l^- 1 

hold for every iVeN and every z satisfying > > 5. 

Proof. Again we start with stating some known result about the Plan- 
cherel-Rotach asymptotics for the Hermite polynomials in the complex 
plane (Theorem 7.185 in [5], see also [18] and the references therein). 
In our notation, uniformly in z from every compact set contained in 
(C U {oo}) \ V2], we have 

(3.30) 

where 

v ; \z + 2 J z-2 

From this relation we obtain for the monic orthogonal Hermite poly- 
nomials (Hn)n>o with respect to the weight exp (— x 2 /2) that 

(3-31) 

H N (zVN) = (2Nf/ 2 H N (z/V2, (2N)- 1 ) 



= 2^N^U(z) «P(^- V5^4)V8) f 1 + fl_ 

uniformly for z from every compact set contained in (CU {oo}) \ [—2, 2]. 
Passing to the normalized polynomials we see that 

(3.32) 



PN (zVN) = (2^1/4(^1)1/2 ^^) 



2 N N N / 2 U(z) exp (N(z - V^M) 2 /8) / f 1_ 
2(27r) 1 /4(AT!)i/2 (7T7p=f^ ^ + \N 



and, in view of the Stirling formula N\ = y/2irN N+ 2e- N (l + 0(1/N)), 
we obtain 

(3.33) 

^ 1 ; 2(27r) 1 /2(AT)i/4( z _yi2-34)7v ^ ^ 
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so that we have 
(3.34) 



2(2ir)'/2(JV)i/4(( z _ V3^4)/2) N V 1^ 

uniformly in z from any compact subset of (C U {oo}) \ [—2, 2]. 

For every complex number z G C \ [—2, 2] let x = x(z) be the root of 

the equation x + x' 1 = z satisfying |x| < 1. Then we have 

x = x(z) = (z - V-2 2 - 4) /2 

where z i— > /c(z) = V-? 2 4 is, by definition, a univalent analytic func- 
tion in C\ [— 2, 2] satisfying the relations k 2 (z) = z 2 — 4 and > for 
t > 2. For every r G (0, 1) the equation \x(z)\ = r and the inequality 
\x(z)\ < r define the ellipse 

(3.35) (^zfKr- 1 + r) 2 + (Qz) 2 /^- 1 - rf = 1 

with focal points —2, 2 and the closed exterior domain E r of this ellipse, 
respectively. Thus for z G E r we have 



|x(z)| = | (2- v^ 2 -4)/2| < r < 1, 

and 

(3.36) 3?((^-v / i 2 ^4)/2) 2 <r 2 <l. 
Note that for every r G (0, 1) 

min \z\ = r _1 — r 

\x(z)\=r 

which for z G E r implies 

2 1 



(3.37) 



-^fz r ^l\ \x(z)\ 



\z - x{z)\ < \z\ + 1 < (1 + 



< 1 + — ; )\z\. 



Further, for z G E r (0 < r < 1) we can write 

\z-2\ = \(x(z) - l) + (x^iz) - l)| = 2\x(z) - - 1| 

= 2|rc(^)|- 1 |l - x(z)\ 2 > 2r- 1 |l - x(z)\ 2 > 2r" 1 |l - rf 

= 2(r -l/2_ r l/2 )2) 

and analogously 

|z + 2| >2(r- 1 / 2 -r 1 / 2 ) 2 . 
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This implies that the function 

\U(z)\ = 



z-2\ 1/A /^ + 2^ 1/4 



z + 2 



is bounded above on E r by a constant depending on r G (0, 1). Applying 
this estimate along with (I3.36P and (I3.37P to the relation (I3.34p we 
arrive, for a fixed r G (0, 1) and every z G E r , at 

(3.38) \pn(zVn)\< C^N-y^ir^zf. 

Now we set z n = Z\J N/ (N — 1) and note that, by convexity of C \ E r , 
zn G E r if z G -E r . Hence, applying (13.381) to Pn-i, we see that for 
every z G E r 

\p N -i{zVN)\ = ^(znVN^T)] < dir^N -lrVV-^lzNl"- 1 

= d(r)(l + l/(iV - l^^N-^L"- 1 ^"- 1 . 
Since (l + 1/(JV - i)) (2Af ~ 3)/4 _ e 1 / 2 as N -> oo, we conclude that 

(3.39) Ipjv-i^VaOI < CaWCJV)- 1 ^- 1 ^)^^- 1 . 
Observe that for every r G (0, 1) the inequality 

Qz > r _1 — r 
implies, by equation (I3.35p . that z G E r . Finally set 



6 /<5 2 n 

r = -2 + VT + 1 

so that r is the solution of the equation r _1 — r = 5 satisfying < r < 1. 
We can now conclude from (I3.38P and (13.391) that for every 5 > with 
some constants C(5) and M(S) the inequalities 

(3.40) |pivOzViV)|< C(5)N- 1 ^M N (5)\z\ N 
and 

(3.41) Ipjv-i^ViV)! < C{5){N)- 1/i M N -\5)\z\ N - 1 

hold for every natural N and every z with \^sz\ > 5. We pass now 
from Hermite polynomials to Hermite functions. Note that for every z 
satisfying \dtz\ > we have 

(3.42) | exp (-z 2 /A)\ < exp ((-3? 2 ^ + % 2 z)/4) < 1 
and, therefore, 

(3-43) \<p k (z)\ < \p k (z)\ 

for every fceE Applying (I3.40p and (13.411) . the Lemma follows. □ 
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Corollary 2. For every real numbers 5 > and A > there exist such 
constants C(5,A) and M(5) that the inequalities 

(3.44) \(p N ((u + t/N)d(H)VN)\< C{5, A)M N \5)\H\ N 
and 

(3.45) \ip N ^((u + t/N) d(H)VN)\< C(5,A)M N ~ 1 (5)\H\ N - 1 

hold for every numbers u £ [—2,2], t £ [—A, A], iV £ N and H £ M 
swc/i t/iat |9f((u + t/N)d(H)) \ > S. 

Proof. Set w = (« + £/iV) d(if). Since 

md(H))\ = y/ (VT+IP - l)/2 < |if|/2, 

we have that 

< (2 + (A/JV))|9f(d(iO)l < (1 + (A/2A^))|iJ| 

and 

\%v)\ N < (1 + (A^JV))^!^^ < C^tf^, 

since 

(1 + (i4/2JV))" - exp(A/2). 
Thus, we proved that 

(3.46) \Q(v)\ N <C(A)\H\ N . 

In view of A3. 18j) . \3lv\ > \$sv\, and we can apply Lemma [5j Omitting 
jy- 1 / 4 j n the bound, the proof is completed. 

□ 

Now we will derive upper bounds for {N^^-Kn^, -, TV^ 1 ) | valid for 
various combinations of the values of its arguments. 

Lemma 6. For every real a £ (0, 1) and A > there exists a constant 
Ri(a,A) > such that for every real u £ [—2 + a, 2 — a], ti,t 2 £ 
[—A, A], H, and every N £ N we have 
(3.47) 

|iV- 1 ^ iV ((w + *i/iV)d(iy), (w + * 2 /iV)d(iy), A^- 1 )! < (^(a, A)) 2JV 
whenev er 

(3.48) |9f(« + U/N)d{H)\ < (3(a/4) 

holds for i = 1,2. 
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Proof. Set 

Ui =u + U/N, Vi = (u + U/N)d(H) for i = l, 2 

and 

u—(u\ + u 2 )/2, v = (vi + v 2 )/2. 

Let the number N = N (a,A) be large enough such that the in- 
equality 

(3.49) A/N < a/2 

holds for every N > N . Observe that by (ET49]) m G [-2 + a/2,2 - 
a/2], i = 1, 2, for N > N . Then it follows from (l3~i8l) and Corollary □ 
that Vi,V2 G Sa/^pfa/q. Therefore, t> G S a /^p^ a /^, too. 
First assume now that we have 

(3.50) \vt -v 2 \ < 1/N. 

Then Remark [1] (with A = 1) applies, and it is clear that, for v\, v 2 , v = 
(t>i + t^)/2 G S a uM a /4\ with fi,f2 satisfying (13.501) . the function 
| sin(A^7r(t> 1 — v 2 )w(vj) / (Nitivi — v 2 )j\ is bounded above by a con- 
stant depending on a and A Therefore, by Remark [TJ this applies to 
the left hand side of ( 13. 47ft as well. 

Now, instead of (I3.50p . we assume that the inequality \v\ — v 2 \ > 1/N 
holds. Then, by the Christoffel-Darboux formula (13.91) . we have 

\N- 1 K N (v 1 ,v 2 ,N- 1 )\ = N 1/2 \K n (Vn Vi ,VNv 2 )\ 

< \Zn(\<pn{VNv 1 )\\<p N -. 1 {VNv 2 )\ + \<Pn{VNv 2 )\\(Pn-i{VNv 1 )\). 

Applying Lemma H] and enlarging, if necessary, the constant in the 
bound, the Lemma is proved. 

□ 

Lemma 7. For every a G (0, 1) and A > there exists a constant 
R 2 (a,A) such that for every u G [—2 + a, 2 — a], ti,t 2 G [—A, A], 
H 6R and N G N we have the relation 
(3.51) 

iN^KNdu+ti/NW^^u+h/NWH)^- 1 )] < T(N)(R 2 (a,A)) N H" 

whenever none of the following pairs of inequalities 

\mu + U/N)d(H))\ </?(a/4), z = l,2, 
\m U + ti/N)d{H))\ >/3(a/4)/2, 2 = 1,2, 
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Proof. Set Ui = (u + t{/N),Vi = Uid(H), i — 1, 2. Since (13.521) does not 
hold we have 1^(^)1 < /?(a/4)/2 and 1 9^(^)1 > P( a /^) f° r some choice 
of i, j G {1, 2}, i 7^ j. This implies 

|u 2 - vil >P{a/4)/2, 
and by the Christoffel-Darboux formula (I3.9P we have 

il^^^^iV- 1 )! = ^1^(^1,^2)1 

^ a/^ A d \'PN(VNv 1 )\\ip N - 1 (VNv 2 )\ + \(p N {VNv 2 )\\(p N ^ 1 {\/Nv 1 )\). 
P(a/4) 

Assume now that, for example, ^(fi)! < /3(a/4)/2 and |S(t>2)| > 
/3(a/4). In addition we have that u x G [-2 + a/2, 2 - a/2] for iV > 
N(a, A). Then, by assertion 2 of Lemma[21 vi G 5 , Q /4 i/ g( Q ,/ 4 )/2 C S , Q ,/ 4j/ g( Q ,/4) , 
and we may apply Lemma i (with the parameter value a/4 instead 
of a) to bound <Pn(V~Nvi) and (fiN-i(VNvi). Furthermore, we apply 
Corollary [2] with 5 = (3{a/A) to bound <Pn(VNv2) and ¥>n-i(VNv2). 
Therefore, we arrive at the inequality 

^\K N ((u + tx/N)d(H), (u + t 2 /N)d(H), l/N) \ 



N 



< C(a, A)V NT{N)M N {f3{a/ 4))H N . 



Replacing the quantity M(/3(a/4)) by an appropriate larger constant, 
we reduce this inequality to the form (I3.5ip . Moreover, this way the 
bound can be made valid for all N G N, which proves the Lemma. 

□ 

Lemma 8. For every 5 > and certain constants C(5),M(5) the 
inequality 

(3.53) l|^(z 1 ,z 2 ,i)|<C'(5)r(Ar)it!^(5)(^ 1 )r(^ 2 ))^ 
holds for every complex numbers Zi,z 2 such that \^sZi\ > 5 and > 

(i = 1, 2), and every natural number N > 5~ 2 . 

Proof. First we consider the case when %tz\$tz2 > 0. Let 

fa >P (x) = ^T«V _1 e" ra P 0,x > 0, 

r(p) 

be the gamma density [7] with the parameters a > 0, p > 0, and let 
P > and iV > 1 be some integers. Then for every complex number z 
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satisfying Jfc > \$$z\ > we have 



(3.54) 

poo pOO 

/ \ z + 9 ^P e -NSt(z+0)V2 d9 / |£ + iQ « z |2P e -S R(€+i9*) a d £ 
/*oo poo 
J\SSz\ J\$tz\ 

POO 

M^r / (e 2 + ( ^ )2) p e -f (? 2 +( ^) 2 )^ 



iV p+1 |3z| 



^ + (^) 2 )/2) Vf ^)^(iV(e 2 + (^) 2 )/2) 

/ ^e-^r/ = T (P + 1) / /y.fa)*? 



iv^i^r^ 1 ^~ 1 i! + "' + p! 

2 p r(p + i)^ iv(^) 2 | [iv(^) 2 ] p 



where we used a well known formula ([7], p. 11) while integrating fi p. 
If |9fz| > 5 then N{^sz) 2 > 1 for iV > 5" 2 , and the bound fl3T54l yields 



(3.55) / \z + 9\ 2P e- N ^ z+ ^/ 2 d9 < e2 p N- 1 T{P + l)|3f( 



oo 



whenever > \^sz\ > 0. In particular, for P = N and P = N — 1 we 
have 



poo 

(3.56) / \z + e\ 2N e- N ^ z+e)2/2 d6 < e2 N N~ l T(N + 
Jo 



zM 2 ^ 1 



and 



POO 

(3.57) / \z + 9\ 2 ^e~ N ^ z+e ^/ 2 d9 < e2 N ~ 1 N- 1 T(N)\Q(z)\ 
Jo 



2N-3 
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Thus, in view of representation (I3.12p and the bounds in Lemma for 
Zi, z 2 such that > 6 and Sftzj > \$$Zi\ (i = 1, 2) we have 

(3.58) 

/■oo 

= (2iV)- 1 / 2 / (< Pn (VNz 1 + t)< Pn „ 1 (^z 2 + t) 
Jo 

+ ifN-iiVNz! + T)^ N {\fNz 2 + r))dr 

POO 

= 2-^ / (i PN (VN(z 1 + e))i fN ^(VN(z 2 + e)) 

Jo 

+ ip N -x(</N(z x + 6))ip N {VN{z 2 + 9)))d6 

POO 



z 2 + 9\ iy e 



+ \ Zl + gf-i e -m*+0) a /4\ 
<C 2 (6)N~ 1 / 2 M 2N -\6)x 

(POO POO 
J \ Zl + e\ 2N e- N ^ + ^' 2 d6 J \ Z2 + e\ 2 ^e 



N -Nft(z 2 +e) 2 /A 



)d6 



N$l(z 2 +8) 2 /2 de 



+[ \ Zl + e\ 2 ( N -i) e -mzi+or/2 de / (\ Z2 + e \ 2N e - N ^ z ^ 2 ' 2 de 



1/2 



1/2-, 



Combining this bound with relations (13.561) and (13.571) we find that 
(3.59) 

\K N {z l ,z 2 ,N- 1 )\ e2 1 / 2 C 2 (5)M 2N - 1 (5)2 N T(N)\$(z 1 )\ N \%z 2 )\ N 
N ~ S 2 N 3 / 2 

The bound (13.591) applies as well when Z\,z 2 satisfy \$sZi\ > 5 and 
$tzi < — 1^1 {i — 1)2). Indeed, this follows from (13.111) . 
Finally, we consider Zi,z 2 such that > \$sz\ > 8 (i = 1,2) and 

dtzi $lz 2 < 0. In this case we have 



< 



< 



1 



\z 2 - Zl \ ~ m Zl )\ + \M(z 2 )\ - |3f(^)| + |3f(z2)| 
and, by (M . 

il^fz!,^,^- 1 )! = -^1^(^1,^2)1 



< (25)' 1 



< -^(1^(^011^-1(^2)1 + 1^(^2)11^-1(^1)1). 
ov iV 
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Applying Lemma Owe obtain 
(3.60) 

\N- 1 K N (z 1 ,z 2 ,N- 1 )\ < C(5)N~ 1 M 2N ~ 1 (5)2 N T(N) (%( Zl )) N (%(z 2 )f . 

Setting R(S) = M (5)^/2 and comparing 13.591 and 13.601 with the factor 
N^ 1 omitted, the proof is completed. 

□ 

Corollary 3. Let a G (0, 1) and A be positive numbers. Then there 
exist positive constants C(a, A) and /^(a) such that for every N G N 
and real numbers u G [—2,2], ti, t 2 G [—A, A], and H we have 
(3.61) 

' <T(N){R 3 (a,A)) 2N H 2N 



whenever the inequalities \$s((u + U/N)d(H)) | > (3(a/4)/2, % — 1,2, 
hold. 

Proof. Set Vi = (u + U/N) d(H) for % = 1,2. From (13.181) we conclude 
that IStoil > |9fVi|, i = 1,2. It follows from (KWi that 

\$S(vi)\ N <C(A)\H\ N ,i = 1,2. 

Then by Lemma [8] with 5 = /5(a/4)/2 we have 

l-K N ((u+%)d(H),(u+^)d(H)^) 

< C(/3(a/4)/2, A)T(N)R 2N ((3(a/4)/2)H 2N 

for every N > (/3(a/4)/2) 2 . Substituting i? in the latter inequality by 
a proper larger constant depending also on A, we obtain (13.611) which 
is valid for all N. 

□ 

Proposition 1. For every a G (0, 1) and A > there exists a constant 
R(a, A) such that for every u G [—2 + a, 2 — a], ti, . . . , t n G [—A, A] 
and H el we have the relations 
(3.62) 



N 
and 



±K N (( u + t±)d(H),(u + %)d(H),±) 



<T(N)(R(a,A)(l + H 2 )) 



N 



GUE, 7 



(3.63) jyn 

<n!(r(iV)) n (^(a,A)) niV (i + iJ 2 r( Af+ i). 
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Proof. For every a G (0, 1) and A > 0, every combination of t\ G 
[—A, A], i 2 G [—A, A], _ff G JR., and JVsN satisfies the assumptions of 
at least one of the assertions in Lemma El Lemma [7] and Corollary [3j 
Hence, at least one of these upper bounds is valid. Writing R(a, A) = 
max( (Ri(a, A)) , R 2 (a, A), (R^at, A)) ), we obtain 

< max(( J R 1 (a, A)f N , T(iV) (i? 2 (a, A)) N H N , T(N) (R 3 (a, A)) 2N H 2N ) 

< (R(a,A)) N max(l,T(N)H N ,T(N)H 2N ) 

< (R(a, A)) N msx(l, T(iV)(l + H 2 ) N/ \ r(JV)(l + #T) 
= r(iV)(i?(a,A)) Af (l + iJT, 

and the bound (I3.62p is proved. 

Now we will derive a bound for the correlation function R^ 1 ^' 1 ^^ 1,11 ^^. 
In view of formulas (13.71) and (13. 8p we have 



(3.64) B^'(x u ...,*„) = ^/ 2 <™'>- 1/2 *i, • • • , ^~ V V) 
It follows from (12.51) that 

/o ar\ D GUE,crs/ x / D HSE,u/iV 2 / x , 

(3.65) R n jy {x 1 ,...,x n ) = R nN 'y N 2^ s {u)du 1 

Jo 

where 
(3.66) 



(2™/ 2 s -V2r( m / 2 ))-V™/2)-i exp ( — uj (2s)), if w > 0, 
0,ifu < 0. 



Since (xi, . . . , x n ) i > R n ^ ' (x\ } . . . , x n ) is an entire function, formula 
(I3.64p can be used to analytically continue the function 
R nI ^' as (xi, . . . ,x n ) in the parameter a to the domain 9far > 0. On 
the other hand, an analytic continuation can be obtained from formu- 
las (I3.65P and (I3.66P (we omit the check of this fact requiring routine 
bounds for 7 miS in the parameter s G C). Since these continuations 
obviously agree, both (13.641) and (13.651) are valid for complex s with 

> 0. In particular, we have the identity 
(3.67) 

R GVE, mW H )N){xu _ ^ = d n (H)R 0^m (xid(H)] . . . jXnd(H )). 
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Along with the determinantal representation (13. 7p (which also admits 
an analytic continuation) and inequality (I3.62p . this leads to 

(3.68) 

AT -n\ D GUE,l/((l+iff)iV) / , /», | . / AT\\ 

N \R UiN '(u + t 1 /N,...,u + t n /N)\ 

= N~ n \d(H) \ n Rn U N' 1/N) ( {u + h/N)d(H), ...,(u + h/N)d(H))\ 

< n!|d(#)| n maxo<Ki<n {N^K^u + U/N)d(H), (u + t 3 /N)d(H))\ n 

< n\\d(H)\ n (T(N) (R(a, A)f (1 + H 2 ) N ) n 
= n\(T(N)) n (R(a, A)) nN (l + H 2 ) n ^ N+1 ^ 

since \d(H)\ = (1 + tf 2 ) 1 / 4 by f[3~m 

□ 

Proof of Theorem. As explained in the end of Section 2 we need to 
show that 

(3.69) 



(Nw(u)) n J_Z v y n ' JV V Nw(u)' ' Nw(u) 

1 , /sin7r(tj — ti 
det v J/ 



uniformly in u,t\, . . . ,t n subject to the conditions formulated in the 
the statement of the Theorem. Throughout the proof we will assume 
that these conditions are satisfied and omit the arguments of -R^ E ' 
whenever that possible. The existence of the integral in the left hand 
side of relation (I3.69P will be a consequence of our estimates. 
Note that R nj N ,s {-Ui, -u ri ) = R^'^m, ...,u n ) since 
Kn(— ui, — Uz) = Kn(ui,v,2) by (13.111) . In view of this property it 
suffices to prove (I3.69|) for u > only. Observe that for real h 



<f) N 2(-h) = (j) N 2(h) 

and 



R GUE,l/((l+iW2/W)jV)/ v _ „GUE,l/((l-iW2/JV)JV) 

- rL n,AT v 1' ' ' ' ' u ™/ ~~ - rL n,Af \ u li ■ ■ ■ i u n) 

for u±, . . . , u n 6 R which shows that (I3.69|) is established if we prove 



(Nw(u)) n Jn 

(3.70) V V " ° 



1 , / sin7r(tj — tj) 
det - v 
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It follows from the central limit theorem for densities that 

f°° 1 
(3.71) / $t(<f> N 2{h))dh 



n->oo 2V2tt 

Here 4>n 2 {~) is a prelimiting characteristic function and l/2\^2n is half 
the value of the limiting standard Gaussian density at 0. Thus, to prove 
(I3.70p it suffices to establish the following relation 



f°° 3ft (4> N 2 (h)R < ^' 1/ ^ 1+ihy ^ /N)N) )dh 
Jo 



(Nw(u))» Jo 

( 3 - 72 ) -detf sm f--y y r*(Mh))dh 

=: I(u,N)-J(N) — ► 

N— *oo 

under the same uniformity constraints as above. Omitting the argu- 
ments ti, . . . , t n , we set 



S n (u, H) = det ( 

and 



/sin(7r(t; - tj)d(H)w(ud(H))/w(u)) 
V 7r(t< - tj)d{H) 

/ sin 7r(ti — tj 



S'n = S^O) = det 



For a given a e (0, 1), let H(a) > be a number whose existence is 
guaranteed by Lemma [31 Let e > be an arbitrary positive number, 
and H = H (a, A, e) G (0, H(a)) be small enough to ensure that 



(3.73) \S n (u, H) - S n \ < e2V2/7r 

for all u E [—2 + a, 2 — a], \H\ < H and ti, . . . ,t n satisfying \ti\ < 
A, . . . , \t n \ < A. Such a number exists because of continuity of the sin- 
kernel and continuity near of the function d(-). Performing the change 
of variable H = hy/2/N, we may write 

(3.74) 

1 /"OO 



iV 



poo 

/ ^<p N ,(NH/V2)RZ EAm+lH)N) )dH 
Jo 



V2(Nw(u)) n Jo 

-- —= / ... dH+ — 

V2{Nw{u)) n Jo V2(Nw{u)) n 

=:/i(e,u,JV)+/ 2 (e,u,JV) 



... dH+ / ... dH 

Ho 
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Quite analogously, we have 
(3.75) 

J(N)=S n U((j) N *{h))dh= —7^ n(<f> N *(NH/V2))dH 
Jo v2 Jo 

rHn at a poo 

= :J 1 {e,N) + J 2 {e,N). 
With this notation we need to show that 



" 1 ^{cj) N 2{NH/V2))dH + —^ / U((f) N 2(NH/V2))dH 

v2 Jhq 



(3.76) I x {e,u,N)- J^N) -— 0, 

N— *oo 

(3.77) I 2 M,iV) -— 0, 
and 

(3.78) J 2 {e,N) — > 

uniformly in ?i G [—2 + a, 2 — a] and G [— ^4, A], i — 1, . . . , n. 

By Proposition [T], for every e, we obtain 

(3.79) 
|/ 2 (e,M,iV)| 

^ -n!(r(iv)) n (i?( a ,A)) nJV r(i + HY N+im \<pN>(^r)\dH. 



< 



71 



2{l + L 2 ) K ~ l 



V2(w(u)y 
Since 

\cf )N 2(NH/V2)\ = {l + H 2 )-^ 2 ^ 

and 
(3.80) 

(1 + tf 2 )-^tf < (1 + L2)A ^ jf (1 + H 2 )dH 
for every K > 1, we have, in particular, that 

\<j> N2 (NH/V2)\dH < — — — ^ ?i - T . 

Note that iV 2 /4-™iV-l/4 > 1 whenever N > N = 2n+^jAn 2 + 5/4. 
For every N > No we obtain from (l3~80l with A" = iV 2 /4 - niV - 1/4 
and L = H that 
(3.82) 

I 2 (e,u,N) < ^ n!(r(iV)) n (#(«, A)) n7V * nr - 

l ' ; " 2y/2(w(u)) n yK )} V V ^ (l + ^o 2 ) /4 " nJV " 5/4 
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Since 
(3.83) 



inf 

ue[-2+e,2-el 



and, by the Stirling formula, (1 + Hq) n2 ^ tends to infinity more rapidly 
than any fixed power of T(N) (or, moreover, the N— th power of any 
positive number), we conclude from ( I3.82p that 

(3.84) h(e,u,N) -— 

uniformly in u G [—2 + a, 2 — a] and ti G [—A, A], i = 1, . . . , n. 

The check of (13.781) is even more straightforward. Indeed, using (13.801) 

with K = N 2 /4 and L = H , we have 



l^(e,iV)| 
ATS 



V2 



Ho 



M((j> N 2(NH/V2))dH 



< 



NS n 
V2 



U N 2(NH/V2)\dH 



Ho 



irNSr, 



h 2y/2(l + L 2 )^ 2 /4-i 



y/2 

and (gZBD follows. 

Now we complete the proof by establishing (I3.76p . Recall that we 
omit the arguments of -R^jy^ 1 ■ All bounds and limit transitions are 
uniform with respect to u G [—2 + a, 2 — a] and ti, t2, ■ ■ ■ , t n satisfying 
\t x \ < A,. . . , \t n \ < A. We have 

(3.85) 

\h(e,u,N)- Ji(e,JV)| 



< 



< 



AT 


rH 







A^ 


/•Ho 


\/2 j 





+ 


N 

7! 




(e,tt 



GUE,l/((l+ii?)JV) 
n,N 



(Nw(u)) n 

GVE,l/((l+iH)N) 
n,N 



S n )(f> N 2{NH/V2) 



dH 



R 



(Nw(u)) r 



S n {u,H))<j> N 2{NH/V2) 



dH 



Ho 



{S n (u,H)-S n )<t> N 2{NH/V2) 
n(2)/. „. Art 



dH 



It follows from the determinantal formula (13.71) and Lemma [3] that 



(3.86) 



sup 

0<H<H 



R 



GUE,l/{(l+iH)N) 
n,N 

(Nw(u)) n 



S n {u,H) 



N^oo 



uniformly in u G [—2 + a, 2 — a], and < A (i — 1, . . . , n). 
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Therefore, using (13.811) . we obtain 



< sup 

0<H<H 



R 



GVE,l/((l+iH)N) 
n,N 



(Nw(u)) n 

\4>^2{h)\dh 



- S n (u,H) 



\4> N 2(h)\dh — ► 



1 



since 
(3.87) 

Jo 1 ' N ~>™ 2V2?r 
Finally, we see from (13. 73f) and (I3.87f) that 

Df\e,u,N) < e. 

Since e is arbitrary small, this completes the proof. □ 
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